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Abstract. Recent developments in the study of primordial black holes (PBHs) will be 
reviewed, with particular emphasis on their formation and evaporation. PBHs could 
provide a unique probe of the early Universe, gravitational collapse, high energy physics 
and quantum gravity. Indeed their study may place interesting constraints on the 
physics relevant to these areas even if they never formed. 



1 Introduction 

Hawking's discovery in 1974 that black holes emit thermal radiation due to 
quantum effects was surely one of the most important results in 20th century 
physics. This is because it unified three previously disparate areas of physics 
- quantum theory, general relativity and thermodynamics - and like all such 
unifying ideas it has led to profound insights. Although not strictly an application 
of quantum gravity theory, the theme of this meeting, it might be regarded as 
a conceptual first step in that direction. Also there is a natural link in that the 
final stage of black hole evaporation, when the black hole is close to the Planck 
mass, can only be understood with a proper theory of quantum gravity. 

In practice, only "primordial black holes" which formed in the early Universe 
could be small enough for Hawking radiation to be important. Such a black hole 
will be referred to by the acronym "PBH" , although this should not be confused 
with the acronym for "Physikcentrum Bad Honnef ' , the institute hosting this 
meeting! Interest in PBHs goes back nearly 35 years and some of the history 
of the subject will be reviewed in Section 2. As will be seen, interest was much 
intensified as a result of Hawking's discovery. Indeed, although it is still not 
definite that PBHs ever formed, it was only through thinking about them that 
Hawking was led to his remarkable insight. Thus the discovery illustrates that 
studying something may be useful even if it does not exist! 

Of course, the subject is much more interesting if PBHs did form and their 
discovery would provide a unique probe of at least four areas of physics: the early 
Universe; gravitational collapse; high energy physics; and quantum gravity. The 
first topic is relevant because studying PBH formation and evaporation can im- 
pose important constraints on primordial inhomogeneities, cosmological phase 
transitions (including inflation) and varying-G models. These topics are covered 
in Sections 3, 4 and 5, respectively. The second topic is discussed in Section 6 
and relates to recent developments in the study of "critical phenomena" and the 
issue of whether PBHs are viable dark matter candidates. The third topic arises 
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because PBH evaporations could contribute to cosmic rays, whose energy dis- 
tribution would then give significant information about the high energy physics 
involved in the final explosive phase of black hole evaporation. This is covered in 
Section 7. The fourth topic arises because it has been suggested that quantum 
gravity effects could appear at TeV scale and this leads to the intriguing possi- 
bility that small black holes could be generated in accelerators experiments or 
cosmic ray events. As discussed in Section 8, this could have striking observa- 
tional consequences. Although such black holes are not technically "primordial" , 
this possibility would have radical implications for PBHs themselves. 

2 Historical Overview 

It was realized many years ago that black holes with a wide range of masses could 
have formed in the early Universe as a result of the great compression associated 
with the Big Bang. A comparison of the cosmological density at a time t after 
the Big Bang with the density associated with a black hole of mass M shows that 
PBHs would have of order the particle horizon mass at their formation epoch: 

PBHs could thus span an enormous mass range: those formed at the Planck 
time (10^^'^s) would have the Planck mass (10~^g), whereas those formed at 1 s 
would be as large as lO'^Af©, comparable to the mass of the holes thought to 
reside in galactic nuclei. By contrast, black holes forming at the present epoch 
could never be smaller than about IMq. 

Zeldovich & Novikov first derived eqn (1) but they were really consid- 
ering "retarded cores" rather than black holes and Hawking was the first 
person to realize that primordial density perturbations might lead to gravita- 
tional collapse on scales above the Planck mass. For a while the existence of 
PBHs seemed unlikely since Zeldovich & Novikov |119j had pointed out that 
they might be expected to grow catastrophically. This is because a simple New- 
tonian argument suggests that, in a radiation-dominated universe, black holes 
much smaller than the horizon cannot grow much at all, whereas those of size 
comparable to the horizon could continue to grow at the same rate as it through- 
out the radiation era. Since we have seen that a PBH must be of order the horizon 
size at formation, this suggests that all PBHs could grow to have a mass of or- 
der lO^^M© (the horizon mass at the end of the radiation era). There are strong 
observational limits on how many such giant holes the Universe could contain, 
so the implication seemed to be that very few PBHs ever existed. 

However, the Zeldovich-Novikov argument was questionable since it neglected 
the cosmological expansion and this would presumably hinder the black hole 
growth. Indeed myself and Hawking were able to disprove the notion that PBHs 
could grow at the same rate as the particle horizon by demonstrating that there 
is no spherically symmetric similarity solution which represents a black hole 
attached to an exact Friedmann model via a sound-wave Since a PBH 
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must therefore soon become much smaller than the horizon, at which stage 
cosmological effects become unimportant, we concluded that PBHs cannot grow 



The realization that small PBHs might exist after all prompted Hawking to 
study their quantum properties. This led to his famous discovery j55| that black 
holes radiate thermally with a temperature 



Only black holes smaller than 10^'^g would have evaporated by the present epoch, 
so eqn (1) implies that this effect could be important only for black holes which 
formed before 10~^^s. 

Despite the conceptual importance of this result, it was bad news for PBH 
enthusiasts. For since PBHs with a mass of lO^^g would be producing photons 
with energy of order 100 MeV at the present epoch, the observational limit 
on the 7-ray background intensity at 100 MeV immediately implied that their 
density could not exceed 10^^ times the critical density |101j . Not only did 
this render PBHs unlikely dark matter candidates, it also implied that there 
was little chance of detecting black hole explosions at the present epoch |1U3| . 
Nevertheless, it was realized that PBH evaporations could still have interesting 
cosmological consequences. In particular, they might generate the microwave 
background il20i or modify the standard cosmological nucleosynthesis scenario 
IHHl or contribute to the cosmic baryon asymmetry P| . PBH evaporations might 
also account for the annihilation-line radiation coming from the Galactic centre 
PU) or the unexpectedly high fraction of antiprotons in cosmic rays [73]. PBH 
explosions occurring in an interstellar magnetic field might also generate radio 
bursts |105| . Even if PBHs had none of these consequences, studying such effects 
leads to strong upper limits on how many of them could ever have formed and 
thereby constrains models of the early Universe. 

Originally it was assumed that PBHs would form from initial inhomogeneities 
but in the 1980s attention switched to several new formation mechanisms. Most 
of the mechanisms were associated with various phase transitions that might 
be expected to occur in the early Universe and there was particular interest in 
whether PBHs could form from the quantum fluctuations associated with the 
many different types of inflationary scenarios. Indeed it soon became clear that 
there are many ways in PBHs serve as a probe of the early Universe and, even 
if they never formed, their non-existence gives interesting information |20|. In 
this sense, they are similar to other "relicts" of the Big Bang, except that they 
derive from much earlier times. 

In the 1990s work on the cosmological consequences of PBH evaporations was 
revitalized as a result of calculations by my PhD student Jane MacGibbon. She 
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realized that the usual assumption that particles are emitted with a black-body 
spectrum as soon as the temperature of the hole exceeds their rest mass is too 
simplistic. If one adopts the conventional view that all particles are composed of 
a small number of fundamental point-like constituents (quarks and leptons), it 
would seem natural to assume that it is these fundamental particles rather than 
the composite ones which are emitted directly once the temperature goes above 
the QCD confinement scale of 250 MeV. One can therefore envisage a black 
hole as emitting relativistic quark and gluon jets which subsequently fragment 
into leptons and hadrons |83l85j and this modifies the cosmological constraints 
considerably [HI] 

Over the last decade PBHs have been assigned various other cosmological 
roles. Some people have speculated that PBH evaporation, rather than proceed- 
ing indefinitely, could cease when the black hole gets down to the Planck mass 
In this case, one could end up with stable Planck mass relics, which would 
provide dark matter candidates |7l25l82j . Although most gamma-ray bursts are 
now known to be at cosmological distances, it has been proposed that some 
of the short period ones could be nearby exploding PBHs | 1(JI28| . Solar mass 
PBHs could form at the quark-hadron phase transition and, since some of these 
should today reside in our Galactic halo, these have been invoked to explain the 
microlensing of stars in the Magellanic Clouds |()4lf)f)lll'^ . 

3 PBHs as a probe of primordial inhomogeneities 

One of the most important reasons for studying PBHs is that it enables one to 
place limits on the spectrum of density fluctuations in the early Universe. This 
is because, if the PBHs form directly from density perturbations, the fraction 
of regions undergoing collapse at any epoch is determined by the root-mean- 
square amplitude e of the fluctuations entering the horizon at that epoch and the 
equation of state p = 7p (0 < 7 < 1). One usually expects a radiation equation 
of state (7 — 1/3) in the early Universe. In order to collapse against the pressure, 
an overdense region must be larger than the Jeans length at maximum expansion 
and this is just ^/j times the horizon size. On the other hand, it cannot be larger 
than the horizon size, else it would form a separate closed universe and not be 
part of our Universe |22| . 

This has two important implications. Firstly, PBHs forming at time t should 
have of order the horizon mass given by eqn (1). Secondly, for a region destined 
to collapse to a PBH, one requires the fractional overdensity at the horizon epoch 
6 to exceed 7. Providing the density fluctuations have a Gaussian distribution 
and are spherically symmetric, one can infer that the fraction of regions of mass 
M which collapse is ^Sl 



f3{M) - €{M) exp 




(4) 



where e(M) is the value of e when the horizon mass is M. The PBHs can have 
an extended mass spectrum only if the fluctuations are scale-invariant (i.e. with 
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e independent of M). In this case, the PBH mass spectrum is given by ^S] 

dn/dM = (a - 2){M/A'Q"'M-^npBHPcr^t (5) 

where « lO^^g is the current lower cut-off in the mass spectrum due to 
evaporations, fipBH is the total density of the PBHs in units of the critical 
density (which itself depends on (3) and the exponent a is determined by the 
equation of state: 

■1 + 37' 



1 + 7 / 

a = 5/2 if one has a radiation equation of state (7=1/3), as expected. This 
means that the density of PBHs larger than M falls off as M^^/^, so most of 
the PBH density is contained in the smallest ones. 

Many scenarios for the cosmological density fluctuations predict that e is at 
least approximately scale- invariant but the sensitive dependence of /9 on e means 
that even tiny deviations from scale-invariance can be important. If e(M) de- 
creases with increasing M, then the spectrum falls off exponentially and most of 
the PBH density is contained in the smallest ones. If e(M) increases with increas- 
ing M , the spectrum rises exponentially and - if PBHs were to form at all - they 
could only do so at large scales. However, the microwave background anisotropics 
would then be larger than observed, so this possibilty can be rejected. 

The current density parameter f2pBH associated with PBHs which form at 
a redshift z or time t is related to /3 by jJi 

npsH^(3nn{i + z)^io^f3i^-j ^loi^l^Y^j (7) 

where f2ii « 10^"* is the density parameter of the microwave background and we 
have used eqn (1). The (1 + z) factor arises because the radiation density scales as 
(l + z)'^, whereas the PBH density scales as (l+z)^. Any limit on fipBH therefore 
places a constraint on (3{M) and the constraints are summarized in Fig.l, which 
is taken from Carr et al. 25. . The constraint for non-evaporating mass ranges 
above 10^''^g comes from requiring Qpbh < 1 but stronger constraints are asso- 
ciated with PBHs smaller than this since they would have evaporated by now 
The strongest one is the 7-ray limit associated with the IQ-'^^g PBHs evapo- 
rating at the present epoch |1U1| . Other ones are associated with the generation 
of entropy and modifications to the cosmological production of light elements 
|98|. The constraints below lO^g are based on the (uncertain) assumption that 
evaporating PBHs leave stable Planck mass relics, in which case these relics are 
required to have less than the critical density (7li25i82) . 

The constraints on /3(M) can be converted into constraints on e(M) using eqn 
(4) and these are shown in Fig. 2. Also shown here are the (non-PBH) constraints 
associated with the spectral distortions in the cosmic microwave background in- 
duced by the dissipation of intermediate scale density perturbations and the 
COBE quadrupole measurement. This shows that one needs the fluctuation am- 
plitude to decrease with increasing scale in order to produce PBHs and the lines 
corresponding to various slopes in the e(Af ) relationship are also shown in Fig. 2. 
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Fig. 1. Constraints on /3(M) 
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Fig. 2. Constraints on e(M) 
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4 PBHs as probe of cosmological phase transitions 

Many phase transitions could occur in the early Universe which lead to PBH 
formation. Some of these mechanisms still require pre-existing density fluctua- 
tions but in others the PBHs form spontaneously even if the Universe starts off 
perfectly smooth. In the latter case, (3{M) depends not on e{M) but on some 
other cosmological parameter. 

4.1 Soft equation of state 

Some phase transitions can lead to the equation of state becoming soft (7 << 1) 
for a while. For example, the pressure may be reduced if the Universe's mass is 
ever channelled into particles which are massive enough to be non-relativistic. 
In such cases, the effect of pressure in stopping collapse is unimportant and the 
probability of PBH formation just depends upon the fraction of regions which 
are sufficiently spherical to undergo collapse; this can be shown to be [20! 

/3 = Ome^^/^. (8) 

The value of /3 is now much less sensitive to e than indicated by eqn (4) and most 
of the PBHs will be smaller than the horizon mass at formation by a factor e'^/^. 
For a given spectrum of primordial fluctuations, this means that there may just 
be a narrow mass range - associated with the period of the soft equation of state 
- in which the PBHs form. In particular, this could happen at the quark-hadron 
phase transition since the pressure may then drop for a while 66 . 

4.2 Collapse of cosmic loops 

In the cosmic string scenario, one expects some strings to self-intersect and form 
cosmic loops. A typical loop will be larger than its Schwarzschild radius by the 
inverse of the factor G/i, where fj, is the mass per unit length. If strings play a 
role in generating large-scale structure, Gfi must be of order 10~^. Hawking [S7j 
showed that there is always a small probability that a cosmic loop will get into a 
configuration in which every dimension lies within its Schwarzschild radius and 
he estimated this to be 

p ^ iG^i)-HG^lx)'^-^ (9) 

where x is the ratio of the loop length to the correlation scale. If one takes x 
to be 3, fipBH > 1 for > 10^^, so he argued that one overproduces PBHs 
in the favoured string scenario. Polnarev & Zemboricz '102' obtained a similar 
result. However, fipBH is very sensitive to x and a slight reduction could still 
give an interesting value |17I41I86] . Note that spectrum (5) still applies since the 
holes are forming with equal probability at every epoch. 
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4.3 Bubble collisions 

Bubbles of broken symmetry might arise at any spontaneously broken symmetry 
epoch and various people, including Hawking, suggested that PBHs could form as 
a result of bubble collisions p2t58.78j . However, this happens only if the bubble 
formation rate per Hubble volume is finely tuned: if it is much larger than the 
Hubble rate, the entire Universe undergoes the phase transition immediately and 
there is not time to form black holes; if it is much less than the Hubble rate, 
the bubbles are very rare and never collide. The holes should have a mass of 
order the horizon mass at the phase transition, so PBHs forming at the GUT 
epoch would have a mass of lO'^g, those forming at the electroweak unification 
epoch would have a mass of lO^^g, and those forming at the QCD (quark-hadron) 
phase transition would have mass of around IMq. Only a phase transition before 
10~^'^s would be relevant in the context of evaporating PBHs. 

4.4 Inflation 

Inflation has two important consequences for PBHs. On the one hand, any PBHs 
formed before the end of inflation will be diluted to a negligible density. Inflation 
thus imposes a lower limit on the PBH mass spectrum: 

M > M,mn = Mpi{Trh/Tpi)-^ (10) 

where Tr// is the reheat temperature and Tpi « 10^^ GeV is the Planck temper- 
ature. The CMB quadrupole measurement implies Tpn ~ lO^^GeV, so Mmin 
certainly exceeds Ig. On the other hand, inflation will itself generate fluctuations 
and these may suffice to produce PBHs after reheating. If the inflaton potential 
is V{(l)), then the horizon-scale fluctuations for a mass-scale M are 

e{M) « [V'/y{Ml.y)]H (11) 

where a prime denotes d/dcj) and the right-hand-side is evaluated for the value 
of (j) when the mass-scale M falls within the horizon. 

In the standard chaotic inflationary scenario, one makes the "slow-roll" and 
"friction-dominated" asumptions: 

^ = {MpiV'/Vf « 1, 77 = Mj>iV"/V « 1. (12) 

Usually the exponent n characterizing the power spectrum of the fluctuations, 
\dk\^ « fc", is very close to but slightly below 1: 

n = 1 + 277 w 1. (13) 

bmce e scales as M^^^")/'^, this means that the fluctuations are slightly increasing 
with scale. The normalization required to explain galaxy formation (e « 10~^) 
would then preclude the formation of PBHs on a smaller scale. If PBH formation 
is to occur, one needs the fluctuations to decrease with increasing mass {n > 1) 
and this is only possible if the scalar held is accelerating sufficiently fast: 



V"/V> {1/2){V'/Vf. 



(14) 
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This condition is certainly satisfied in some scenarios |23| and, if it is, eqn (4) 
implies tliat tfie PBH density will be dominated by the ones forming immedi- 
ately after reheating. Since each value of n corresponds to a straight line in Fig. 2, 
any particular value for the reheat time ti corresponds to an upper limit on n. 
This limit is indicated in Fig. 3, which is taken from Carr et al. j25| apart from 
a correction pointed out by Green & Liddle |47j . Similar constraints have now 
been obtained by several other people |15I72| . The figure also shows how the con- 
straint on n is strengthened if the reheating at the end of inflation is sufficiently 
slow for there to be a dust-like phase 021 ■ PBHs have now been used to place 
constraints on many other sorts of inflationary scenarios - supernatural \W^ , su- 
persymmetric hybrid |4(JI68| . oscillating |1HJ| . preheating 9'34 38 50 and 
running mass |79| - as well as a scenarios in which the inflaton serves as the dark 
matter [ST] . 



Bullock & Primack ^B] and Ivanov 63 have questioned whether the Gaus- 
sian assumption which underlies eqn (4) is valid in the context of inflation. So 
long as the fluctuations are small {S(p/(j) << 1), as certainly applies on a galac- 
tic scale, this assumption is valid. However, for PBH formation one requires 
(50/(/) ~ 1, and, in this case, the coupling of different Fourier modes destroys the 
Gaussianity. Their analysis suggests that j3{M) is much less than indicated by 
eqn (4) but it still depends very sensitively on e. 

Not all inflationary scenarios predict that the spectral index should be con- 
stant. Hodges & Blumenthal 'ST pointed out that one can get any form for the 
fluctuations whatsoever by suitably choosing the form of V{(j)). For example, eqn 
(11) suggests that one can get a spike in the spectrum by flattening the potential 
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Fig. 3. Constraints on spectral index n in terms of reheat time ii 
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over some mass range (since the fluctuation diverges when V' goes to 0). This 
idea was exploited by Ivanov et al. ^64 , who fine-tuned the position of the spike 
so that it corresponds to the microlensing mass-scale. 

5 PBHs as a probe of a varying gravitational constant 

The PBH constraints would be severely modified if the value of the gravitational 
"constant" G was different at early times. The simplest varying-G model is 
Brans-Dicke (BD) theory in which G is associated with a scalar field (f> 
and the deviations from general relativity are specified by a parameter w. A 
variety of astrophysical tests currently require |a;| > 500, which implies that the 
deviations can only ever be small However, there exist generalized scalar- 

tensor theories |11I97I11'2] in which oj is itself a function of (f> and these lead to 
a considerably broader range of variations in G. In particular, it permits to to 
be small at early times (allowing noticeable variations of G then) even if it is 
large today. In the last decade interest in such theories has been revitalized as a 
result of early Universe studies. Extended inflation explicitly requires a model in 
which G varies [7S] and, in higher dimensional Kaluza-Klein-type cosmologies, 
the variation in the sizes of the extra dimensions also naturally leads to this 

mm- 

The behaviour of homogeneous cosmological models in BD theory is well 
understood |^. They are vacuum-dominated at early times but always tend 
towards the general relativistic solution during the radiation-dominated era. This 
means that the full radiation solution can be approximated by joining a BD 
vacuum solution to a general relativistic radiation solution at some time which 
may be regarded as a free parameter of the theory. However, when the matter 
density becomes greater than the radiation density at around lO^y, the equation 
of state becomes dustlike {p = 0) and G begins to vary again. 

The consequences of the cosmological variation of G for PBH evaporation 
depend upon how the value of G near the black hole evolves. Barrow ^ intro- 
duces two possibilities: in scenario A, G everywhere maintains the background 
cosmological value (so cj) is homogeneous); in scenario B, it preserves the value it 
had at the formation epoch near the black hole even though it evolves at large 
distances (so cj) becomes inhomogeneous). On the assumption that a PBH of 
mass M has a temperature and mass-loss rate 

T(xiGM)-\ M(x{GM)-^, (15) 

with G = G{t) in scenario A and G = G{M) in scenario B, Barrow & Carr 
PI calculate how the evaporation constraints summarized in Fig.l are modified 
for a wide range of varying-G models. The question of whether scenario A or 
scenario B is more plausible has been studied in several papers |:jll43lh:>lHh| but 
is still unresolved. 



Primordial Black Holes 11 



6 PBHs as a probe of gravitational collapse 

The criterion for PBH formation given in Section 3 is rather simplistic and not 
based on a detailed calculation. The first numerical studies of PBH formation 
were carried out by Nadezhin et al. |92j . These roughly confirmed the criterion 
^ > 7 for PBH formation, although the PBHs could be be somewhat smaller 
than the horizon. In recent years several groups have carried out more detailed 
hydrodynamical calculations and these have refined the ^ > 7 criterion and hence 
the estimate for (3{M) given by eqn (4). Niemeyer & Jedamzik [21] find that one 
needs S > 0.8 rather than 6 > 0.3 to ensure PBH formation and they also find 
that there is little accretion after PBH formation, as expected theoretically |22| . 
Shibata & Sasaki |108j reach similar conclusions. 

A particularly interesting development has been the application of "critical 
phenomena" to PBH formation. Studies of the collapse of various types of spher- 
ically symmetric matter fields have shown that there is always a critical solution 
which separates those configurations which form a black hole from those which 
disperse to an asymptotically flat state. The configurations are described by 
some index p and, as the critical index pc is approached, the black hole mass is 
found to scale as {p — PcY^ for some exponent ry. This effect was first discovered 
for scalar fields [2^1 but subsequently demonstrated for radiation and then 
more general fluids with equation of state p = jp |75l9()j . 

In all these studies the spacetime was assumed to be asymptotically flat. 
However, Niemeyer & Jedamzik |95| have recently applied the same idea to 
study black hole formation in asymptotically Friedmann models and have found 
similar results. For a variety of initial density perturbation profiles, they find 
that the relationship between the PBH mass and the the horizon-scale density 
perturbation has the form 

M = KMh{5 - ScV (16) 

where Mh is the horizon mass and the constants are in the range 0.34 < 7 < 0.37, 
2A< K < 11.9 and 0.67 < 5c < 0.71 for the various configurations. Since M ^ 
as (5 — > (5c, this suggests that PBHs may be much smaller than the particle hori- 
zon at formation and it also modifies the mass spectrum |45I48I76I117| . However, 
it is clear that a fluid description must break down if they are too small and 
recent calculations by Hawke & Stewart show that black holes can only form 
on scales down to 10~^ of the horizon mass. 

There has also been interest recently in whether PBHs could have formed at 
the quark-hadron phase transition at 10~^s because of a temporary softening of 
the equation of state then. Such PBHs would naturally have the sort of mass 
required to explain the MACHO microlensing results 66 . If the QCD phase 
transition is assumed to be of 1st order, then hydrodynamical calculations show 
that the value of 5 required for PBH formation is indeed reduced below the value 
which pertains in the radiation case [HZ| . This means that PBH formation will be 
strongly enhanced at the QCD epoch, with the mass distribution being peaked 
around the horizon mass. One of the interesting implications of this scenario is 
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the possible existence of a halo population of binary black holes • With a full 
halo of such objects, there could then be 10^ binaries inside 50 kpc and some 
of these could be coalescing due to gravitational radiation losses at the present 
epoch. If the associated gravitational waves were detected, it would provide a 
unique probe of the halo distribution |fi2j . 



7 PBHs as a probe of high energy physics 

We have seen that a black hole of mass M will emit particles like a black-body 
of temperature [nE] 

T»10-(^f)"K.(Ji-)"GeV. (IT) 

This assumes that the hole has no charge or angular momentum. This is a 
reasonable assumption since charge and angular momentum will also be lost 
through quantum emission but on a shorter timescale that the mass |100| . This 
means that it loses mass at a rate 

M = -5 X lO'^%M/g)-'^f{M) g s'^ (18) 

where the factor f{M) depends on the number of particle species which are light 
enough to be emitted by a hole of mass M, so the lifetime is 

t{M) = 6 X 10-^'^f{M)-\M/gf s. (19) 

The factor / is normalized to be 1 for holes larger than 10^^ g and such holes 
are only able to emit "massless" particles like photons, neutrinos and gravitons. 
Holes in the mass range 10^^ g < M < 10^^ g are also able to emit electrons, 
while those in the range 10" g < M < 10^5 g emit muons which subsequently 
decay into electrons and neutrinos. The latter range includes, in particular, the 
critical mass for which r equals the age of the Universe. If the total density 
parameter is 1, this can be shown to M, = 4.4 x where h is the 

Hubble parameter in units of 100 [Mj . 

Once M falls below lO^^g, a black hole can also begin to emit hadrons. 
However, hadrons are composite particles made up of quarks held together by 
gluons. For temperatures exceeding the QCD confinement scale of Aqcd = 250— 
300 GeV, one would therefore expect these fundamental particles to be emitted 
rather than composite particles. Only pions would be light enough to be emitted 
below Aqcd- Since there are 12 quark degrees of freedom per flavour and 16 
gluon degrees of freedom, one would also expect the emission rate (i.e. the value 
of /) to increase dramatically once the QCD temperature is reached. 

The physics of quark and gluon emission from black holes is simplified by a 
number of factors. Firstly, one can show that the separation between successively 
emitted particles is about 20 times their wavelength, which means that short 
range interactions between them can be neglected. Secondly, the condition T > 
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ylqcD implies that their separation is much less than ^q^p ~ 10^^'^ cm (the 
characteristic strong interaction range) and this means that the particles are 
also unaffected by strong interactions. The implication of these three conditions 
is that one can regard the black hole as emitting quark and gluon jets of the 
kind produced in collider events. The jets will decay into hadrons over a distance 
which is always much larger than GM , so gravitational effects can be neglected. 
The hadrons may then decay into astrophysically stable particles through weak 
and electomagnetic decays. 

To find the final spectra of stable particles emitted from a black hole, one 
must convolve the Hawking emission spectrum with the jet fragmentation func- 
tion. This gives the instantaneous emission spectrum shown in Fig. 4 for a T = 
1 GcV black hole (SSj. The direct emission just corresponds to the small bumps 
on the right. All the particle spectra show a peak at 100 MeV due to pion decays; 
the electrons and neutrinos also have peaks at 1 MeV due to neutron decays. In 
order to determine the present day background spectrum of particles generated 
by PBH evaporations, one must first integrate over the lifetime of each hole of 
mass M and then over the PBH mass spectrum In doing this, one must 
allow for the fact that smaller holes will evaporate at an earlier cosmological 
epoch, so the particles they generate will be redshifted in energy by the present 
epoch. 




-3 -2 -1 
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Fig. 4. Instantaneous emission from a 1 GeV black hole 

If the holes arc uniformly distributed throughout the Universe, the back- 
ground spectra should have the form indicated in Fig. 5. All the spectra have 
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rather similar shapes: an E^^ fall-ofF for E > 100 MeV due to the final phases 
of evaporation at the present epoch and an tail for E < 100 MeV due to 
the fragmentation of jets produced at the present and earlier epochs. Note that 
the E^^ tail generally masks any effect associated with the mass spectrum of 
smaller PBHs which evaporated at earlier epochs |19j . 

The situation is more complicated if the PBHs evaporating at the present 
epoch are clustered inside our own Galactic halo (as is most likely). In this case, 
any charged particles emitted after the epoch of galaxy formation (i.e. from PBHs 
only somewhat smaller than M*) will have their flux enhanced relative to the 
photon spectra by a factor ^ which depends upon the halo concentration factor 
and the time for which particles are trapped inside the halo by the Galactic 
magnetic field. This time is rather uncertain and also energy-dependent. At 
100 MeV one has ^ ~ 10^ for electrons or positrons and ^ ~ 10'* for protons and 
antiprotons. MacGibbon & Carr [H3 first used the observed cosmic ray spectra 
to constrain I?pbh but their estimates have recently been updated. 




-4 -3 -2 -1 
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Fig. 5. Spectrum of particles from uniformly distributed PBHs 



7.1 Gamma-rays 

Recent EGRET observations jl09) give a 7-ray background of 

,p / p X -2.10±0.03 

^ = 7.3(±0.7)xlO-»(^) c„-GeV- (20) 
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between 30 MeV and 120 GeV. Carr & MacGibbon 24 showed that this leads 
to an upper hmit 

^^PBH < (5.1± 1.3) X 10"^/^"^ (21) 

which is a refinement of the original Page-Hawking hmit, but the form of the 
spectrum suggests that PBHs do not provide the dominant contribution. If PBHs 
are clustered inside our own Galactic halo, then there should also be a Galactic 7- 
ray background and, since this would be anisotropic, it should be separable from 
the extragalactic background. The ratio of the anisotropic to isotropic intensity 
depends on the Galactic longtitude and latitude, the ratio of the core radius to 
our Galactocentric radius, and the halo flattening. Wright claims that such a 
halo background has been detected |114) . His detailed fit to the EGRET data, 
subtracting various other known components, requires the PBH clustering factor 
to be (2 — 12) X 10^h~^, comparable to that expected. 

7.2 Antiprotons 

Since the ratio of antiprotons to protons in cosmic rays is less than 10^** over 
the energy range 100 MeV — 10 GeV, whereas PBHs should produce them in 
equal numbers, PBHs could only contribute appreciably to the antiprotons |111| . 
It is usually assumed that the observed antiproton cosmic rays are secondary 
particles, produced by spallation of the interstellar medium by primary cosmic 
rays. However, the spectrum of secondary antiprotons should show a steep cut- 
off at kinetic energies below 2 GeV, whereas the spectrum of PBH antiprotons 
should increase with decreasing energy down to 0.2 GeV, so this provides a 
distinct signature [75| . 

MacGibbon & Carr originally calculated the PBH density required to explain 
the interstellar antiproton flux at 1 GeV and found a value somewhat larger than 
the 7-ray limit |H1]. More recent data on the antiproton flux below 0.5 GeV comes 
from the BESS balloon experiment |118| and Maki et al. [HHI have tried to fit 
this data in the PBH scenario. They model the Galaxy as a cylindrical diffusing 
halo of diameter 40 kpc and thickness 4-8 kpc and then using Monte Carlo 
simulations of cosmic ray propagation. A comparison with the data shows no 
positive evidence for PBHs (i.e. there is no tendency for the antiproton fraction 
to tend to 0.5 at low energies) but they require the fraction of the local halo 
density in PBHs to be less than 3 x 10^* and this is stronger than the 7-ray 
background limit. A more recent attempt to fit the observed antiproton spectrum 
with PBH emission comes from Barrau et al. [Hj and is shown in Fig. 6. A key test 
of the PBH hypothesis will arise during the solar minimum period because the 
flux of primary antiprotons should be enhanced then, while that of the secondary 
antiprotons should be little affected |91j . 

7.3 PBH Explosions 

One of the most striking observational consequences of PBH evaporations would 
be their final explosive phase. However, in the standard particle physics picture, 
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PBH antiproton spectra 



in 
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Fig. 6. Comparison of PBH emission and antiproton data from Barrau et al. 



where the number of elementary particle species never exceeds around 100, the 
likelihood of detecting such explosions is very low. Indeed, in this case, observa- 
tions only place an upper limit on the explosion rate of 5 x 10*pc~'^y~^ |1I1U7| . 
This compares to Wright's 7-ray halo limit of 0.3 pc~^y~^ and the Maki et al. 
antiproton limit of 0.02 pc~^y~^. 

However, the physics at the QCD phase transition is still uncertain and the 
prospects of detecting explosions would be improved in less conventional particle 
physics models. For example, in a Hagedorn-type picture, where the number of 
particle species exponentiates at the the quark-hadron temperature, the upper 
Hmit is reduced to 0.05 pc~'^y~^ Chne and colleagues have argued that one 
might expect the formation of a QCD fireball at this temperature and this 
might even explain some of the short period 7-ray bursts observed by BATSE 
|28| . They claim to have found 42 candidates of this kind and the fact that their 
distribution matches the spiral arms suggests that they are Galactic. Although 
this proposal is speculative and has been disputed it has the attraction 
of making testable predictions (eg. the hardness ratio should increase as the 
duration of the burst decreases). A rather different way of producing a 7-ray 
burst is to assume that the outgoing charged particles form a plasma due to 
turbulent magnetic field effects at sufficiently high temperatures |10| . 
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Some people have emphasized the possibility of detecting very high energy 
cosmic rays from PBHs using air shower techniques '3 1'5 1"77' . However, re- 
cently these efforts have been set back by the claim of Heckler that QED 
interactions could produce an optically thick photosphere once the black hole 
temperature exceeds Tcrit — 45 GeV. In this case, the mean photon energy is 
reduced to me{TBH/Tcrity^^ , which is well below Tbh, so the number of high 
energy photons is much reduced. He has proposed that a similar effect may op- 
erate at even lower temperatures due to QCD effects |HD]- Several groups have 
examined the implications of this proposal for PBH emission |29lf)9j . However, 
these arguments should not be regarded as definitive since MacGibbon et al. 
claim that QED and QCD interactions are never important |87| . 

8 PBHs as a probe of quantum gravity 

In the standard Kaluza-Kclin picture, the extra dimensions are assumed to be 
compactified on the scale of the Planck length. This means that the influence of 
these extra dimensions only becomes important at an energy scale of lO^^GeV 
and this is also presumbaly the scale on which quantum gravity effects become 
significant. In particular, such effects are only important for black hole evapo- 
rations once the black hole mass gets down to the Planck mass of 10~^g. Con- 
ceivably, this could result in black hole evaporation ceasing, so that one ends up 
with stable Planck-mass relics, and this leads to the sort of "relics" constraints 
indicated in Fig.l, Fig. 2 and Fig. 3. Various non-quantum-gravitational effects 
(such as higher order corrections to the gravitational Lagrangian or string ef- 
fects) could also lead to stable relics but the relic mass is always close to 
the Planck mass. 

In "brane" versions of Kaluza-Klein theory, some of the extra dimensions can 
be much larger than the Planck length and this means that quantum gravity 
effects may become important at a much smaller energy scale. If the internal 
space has n dimensions and a compact volume Vn, then Newton's constant Gn 
is related to the higher dimensional gravitational constant Gd and the value of 
the modified Planck mass Mp is related to the usual 4-dimensional Planck mass 
M4 by the order-of-magnitude equations: 

Gn - Go/Vn, m;+2 ~ M|/K. (22) 

The same relationship applies if one has an infinite extra dimension but with 
a "warped" geometry, provided one interprets Vn as the "warped volume" . In 
the standard model, Vn ~ 1/M" and so Mp ^ M4. However, with large extra 
dimensions, one has Vn >> l/M'2 and so Mp << M4. In particular, this might 
permit quantum gravitational effects to arise at the experimentally observable 
TeV scale. 

If this were true, it would have profound implications for black hole formation 
and evaporation since black holes could be generated in accelerator experiments, 
such as the Large Hadron Collider (LHC). Two partons with centre-of-mass 
energy ^/s will form a black hole if they come within a distance corresponding 
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to the Schwarzschild radius rs for a black hole whose mass Mbh is equivalent 
to that energy |33l42lldH| . Thus the cross-section for black hole production is 

aBH«vrr|e(V^-Af™") (23) 

where M]^™ is the mass below which the semi-classical approximation fails. Here 
the Schwarzschild radius itself depends upon the number of internal dimensions: 



1 (MbhV'^'^''^ 



so that gbh s1/("+i). This means that the cross-section for black hole produc- 
tion in scattering experiments goes well above the cross-section for the standard 
model above a certain energy scale and in a way which depends on the number 
of extra dimensions. 

The evaporation of the black holes produced in this way will produce a char- 
acteristic signature |33I42I10B| because the temperature and lifetime of the black 
holes depend on the number of internal dimensions: 

Thus the temperature is decreased relative to the standard 4-dimensional case 
and the lifetime is increased. The important qualitative effect is that a large 
fraction of the beam energy is converted into transverse energy, leading to large- 
multiplicity events with many more hard jets and leptons than would otherwise 
be expected. In principle, the formation and evaporation of black holes might 
be observed by LHC by the end of the decade and this might also allow one to 
experimentally probe the number of extra dimensions. On the other hand, this 
would also mean that scattering processes above the Planck scale could not be 
probed directly because they would be hidden behind a black hole event horizon. 

Similar effects could be evident in the interaction between high energy cos- 
mic rays and atmospheric nucleons. Nearly horizontal cosmic ray neutrinos would 
lead to the production of black holes, whose decays could generate deeply pen- 
etrating showers with an electromagnetic component substantially larger than 
that expected with conventional neutrino interactions. Several authors have stud- 
ied this in the context of the Pierre Auger experiment, with event rates in excess 
of one per year being predicted |i2i36al06j . Indeed there is a small window of 
opportunity in which Auger might detect such events before LMC. 

It should be stressed that the black holes produced in these processes should 
not themselves be described as "primordial" since they do not form in the early 
Universe. On the other hand, it is clear that the theories which predict such 
processes will also have profound implications for the formation and evaporation 
of those black holes which do form then. This is because, at sufficiently early 
times, the effects of the extra dimensions must be cosmologically important. 
However, these effects are not yet fully understood. 
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9 Conclusions 

We have seen that PBHs could provide a unique probe of the early Universe, 
gravitational collapse, high energy physics and quantum gravity. In the "early 
Universe" context, particularly useful constraints can be placed on inflationary 
scenarios and on models in which the value of the gravitational "constant" G 
varies with cosmological epoch. In the "gravitational collapse" context, the ex- 
istence of PBHs could provide a unique test of the sort of critical phenomena 
discovered in recent numerical calculations. In the "high energy physics" con- 
text, information may come from observing cosmic rays from evaporating PBHs 
since the constraints on the number of evaporating PBHs imposed by gamma-ray 
background observations do not exclude their making a significant contribution 
to the Galactic flux of electrons, positrons and antiprotons. Evaporating PBHs 
may also be detectable in their final explosive phase as gamma-ray bursts if suit- 
able physics is invoked at the QCD phase transition. In the "quantum gravity" 
context, the formation and evaporation of small black holes could lead to ob- 
servable signatures in cosmic ray events and accelerator experiments, providing 
there are extra dimensions and providing the quantum gravity scale is around a 
TeV. 
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